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ABSTRACT 

We report on the observation of 7-rays above 25GeV from the Crab pulsar (PSR B0532+21) using 
the MAGIC I telescope. Two data sets from observations during the winter period 2007/2008 and 
2008/2009 are used. In order to discuss the spectral shape from 100 MeV to 100 GeV, one year of 
public Fermi Large Area Telescope (Fermi-LAT) data are also analyzed to complement the MAGIC 
data. The extrapolation of the exponential cutoff spectrum determined with the Fermi-hAT data is 
inconsistent with MAGIC measurements, which requires a modification of the standard pulsar emission 
models. In the energy region between 25 and 100 GeV, the emission in the PI phase (from —0.06 to 
0.04, location of the main pulse) and the P2 phase (from 0.32 to 0.43, location of the interpulse) can 
be described by power laws with spectral indices of — 3.1± 1.0 s t a t±0.3 sys t and — 3.5±0.5 s t a t±0.3 S!/S (, 
respectively. Assuming an asymmetric Lorentzian for the pulse shape, the peak positions of the main 
pulse and the interpulse are estimated to be at phases —0.009 ± 0.007 and 0.393 ± 0.009, while the 
full widths at half-maximum are 0.025 ± 0.008 and 0.053 ± 0.015, respectively. 

Subject headings: Gamma rays: stars — pulsars: individual (Crab pulsar, PSR B0531+21) 
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1. INTRODUCTION 

The pulsar B0531+21, also commonly known as the 
Crab pulsar, is the compact object left over after a his- 
toric supernova explosion that occurred in the year 1054 
AD. Its energetic pulsar wind creates a pulsar wind neb- 
ula, the Crab Nebula, an emitter of strong and steady 
radiation. The pulsar and pulsar wind nebula have been 
observed and studied in almost the entire accessible elec- 
tromagnetic spectrum from about 10~ 5 eV (radio emis- 
sion) to nearly 100 TeV (very high energy 7-rays, e.g., 
Aharonian et al. (2004)). The nebular emission is com- 
monly used as a standard candle for astronomy in various 
energy ranges. Recently, 7-ray flares from the Crab Neb- 
ula were discovered in the GeV range (Tavani et al. 2011; 
Abdo et al. 2011) and a hint of increased flux in the TeV 
range during a GeV flare was also reported (ATel 2921). 

The Crab pulsar and several other pulsars are 
amongst the brightest known sources at 1 GeV. How- 
ever, a spectral steepening made their detection above 
10 GeV elusive despite numerous efforts (e.g., Aharonian 
et al. (2007); de Naurois et al. (2002); Lessard et al. 
(2000)). The energy thresholds of imaging atmospheric 
Cherenkov telescopes (IACTs) were, in general, too high, 
while the 7-ray collection area of satellite-borne detectors 
was too small to detect pulsars above 10 GeV. 

On the other hand, a precise measurement of the en- 
ergy spectrum at and above the steepening is an impor- 
tant test for the standard pulsar models, such as the po- 
lar cap (PC), outer gap (OG), and slot gap (SG) models. 

In the PC model, emission takes place within a few 
neutron star (NS) radii above a PC surface (Arons & 
Scharlemann 1979; Daugherty & Harding 1982, 1996). 
There, high energy gamma-rays above ~ 1 GeV should be 
absorbed by a strong magnetic field (Baring 2004), which 
results in a very sharp cutoff (so-called super-exponential 
cutoff) in the energy spectrum. 

Extending the original idea by Arons (1983), Muslimov 
& Harding (2003, 2004b,a) and Dyks et al. (2004) inves- 
tigated the possibility of high-energy emission along the 
flaring field lines at high altitudes. This type of emis- 
sion, a SG emission, can be observed at all viewing angle 
and for most cases emission from the two poles can be 
observed. The SG model predicts an exponential cutoff 
above 1 GeV (Harding et al. 2008). However, such geo- 
metrically thin emission models reproduce less than 20 % 
of the observed 7-ray fluxes of the Crab and Vela pulsars 
(Hirotani 2008). 

Seeking a different possibility of high-altitude emis- 
sions, Cheng et al. (1986a, b) proposed the OG model, 
hypothesizing that the emission zone is located in higher 
altitudes, beyond the so-called null-charge surface. Sub- 
sequently, Romani & Yadigaroglu (1995) and Romani 
(1996) developed the caustic model of the OG emissions. 
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A three-dimensional version of such a geometrical model 
of OG emissions was investigated by Cheng et al. (2000). 
Recently, Romani & Watters (2010) presented an atlas of 
pulse properties and proposed a method to discriminate 
different emission models from the geometrical point of 
view. 

It is noteworthy that all the existing OG or SG models 
predict that the highest-energy photons are emitted via 
curvature radiation and that an exponential cutoff ap- 
pears around 10 GeV in the spectrum of the Crab pulsar 
(e.g., (Tang et al. 2008)). 

The MAGIC I telescope in its standard trigger mode 
has the worldwide lowest threshold of all currently op- 
crating IACTs, around 60 GeV. A previous study of 
MAGIC data above 60 GeV revealed a 2.9 o excess from 
the Crab pulsar (Albert et al. 2008c). Following up on 
this hint, we investigated an alternative trigger concept, 
the sum trigger (see section 2.2), which lowered the en- 
ergy threshold of MAGIC to about 25 GeV. Using this 
new trigger, we observed the Crab pulsar between 2007 
October and 2008 February and detected high-energy 7- 
ray emission from the Crab pulsar with a significance 
of 6.4er (Aliu et al. 2008). This detection suggests the 
distance of the emission region from the stellar surface 
to be larger than 6.2 ±0.2 stat ± 0.3 sys t times the stellar 
radius, which ruled out the PC model as a viable expla- 
nation of the observed emission. This initial detection 
has been briefly reported in Aliu et al. (2008). In win- 
ter 2008/2009, the Crab pulsar was observed again with 
MAGIC using the sum trigger. 

In 2008 August, the new satellite-borne 7-ray detec- 
tor with 1 m 2 collection area, the Fermi Large Area 
Telescope (Fermi-LAT), became operational and mea- 
sured the spectra of 7-ray pulsars up to a few tens of 
GeV (Abdo et al. 2010c). All the energy spectra are 
consistent with a power law with an exponential cutoff, 
though statistical uncertainties above 10 GeV are rather 
large 34 . The cutoff energies are typically between 1 GeV 
and 4 GeV. These Fermi-LAT measurements also disfa- 
vor the PC model and support the OG and the SG model 
(Abdo et al. 2010c). 

However, the cutoff energy of the Crab pulsar deter- 
mined with the Fermi-LAT under the exponential cutoff 
assumption is ~ 6 GeV, an unlikely value for the signal 
above 25 GeV detected by MAGIC. In order to verify 
the exponential cutoff spectrum, a precise comparison of 
the energy spectra measured by the two instruments is 
needed. The recent detection of the Crab pulsar above 
100 GeV by the VERITAS Collaboration has shown that 
indeed the energy spectrum above the break is not con- 
sistent with an exponential cutoff but that it is better 
described by a broken power law (Aliu et al. 2011). It 
is, however, not clear whether the spectrum continues as 
a power law after the break or there is another compo- 
nent above 100 GeV in addition to the exponential cutoff 
spectrum because of missing flux measurements in the 
intermediate energy range from 25 GeV to 100 GeV. 

The main objectives of this paper are the evaluation 
of the exponential cutoff spectrum of the Crab pulsar 

34 For some of the pulsars, the phase averaged spectrum devi- 
ates from the exponential cutoff but phase-resolved analyses re- 
vealed that the spectrum of each small pulse-phase interval is still 
consistent with the exponential cutoff (Abdo et al. 2010d,b). 
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with the MAGIC data and the presentation of its energy 
spectrum between 25 GeV and 100 GeV. We also give de- 
tails of the MAGIC observations, the data selection, the 
analysis and physics results. We report, for the first time, 
separate energy spectra and a pulse profile analysis above 
25 GeV for both the main pulse and the interpulse. The 
large majority of the results presented in this paper are 
extracted from the PhD thesis of Takayuki Saito (Saito 
2010). The paper has the following structure: After de- 
scribing the MAGIC telescope and the sum trigger in 
Section 2, we present the observation and details of the 
data processing in Section 3. The detection of pulsed 
emission is described in Section 4. Based on the MAGIC 
detection and the Ferraz-LAT measurements, the evalu- 
ation of the exponential cutoff assumption is reported in 
Section 5. Energy spectra for the main pulse, the inter- 
pulse, and the summed pulsed emission are presented in 
Section 6, followed by a discussion of the pulse profile in 
Section 7. We conclude the paper in Section 8 together 
with a theoretical interpretation of the spectrum and an 
outlook on what can be expected in the near future. 

2. THE MAGIC I TELESCOPE 
2.1. System Overview 

The MAGIC I telescope is a new generation IACT lo- 
cated on the Canary island of La Palma (27?8 N, 17?8 
W, 2225 m asl). Part of the Chercnkov photons emitted 
by charged particles in an air shower are collected by a 
parabolic reflector with a diameter of 17 m and focused 
onto a fine pixclizcd camera, providing an image of the 
air shower. The camera comprises 577 photomultiplicr 
(PMTs) and has a field of view of ~ 3?6 diameter. 

The fast analog PMT signals are transported via op- 
tical fibers to the counting house where the signals are 
processed by the trigger system and recorded by the data 
acquisition (DAQ) system. The trigger is normally de- 
rived from the pixels in the innermost camera area (the 
trigger area) of around 1° radius (325 pixels). Each sig- 
nal from the pixels in the trigger area is amplified and 
split into two signals equal in amplitude. The signals are 
routed to the trigger system and to the DAQ system, re- 
spectively. Signals from the non-trigger pixels enter the 
DAQ system directly. The trigger criteria of the standard 
digital trigger (hereafter called standard trigger) system 
are applied in two steps: each optical signal from the 
trigger area is converted to an electric one and examined 
by a discriminator with a computer-controlled threshold 
level; the threshold level is typically 6 — 7 photoelectrons 
(PhEs). The digital signals are then processed by a topo- 
logical pattern logic, which searches for a close-packed 
group of four compact-next-neighbor pixels firing within 
a time window of ~ 6 ns. 

In the standard trigger concept, only signals above the 
preset threshold in four compact-ncxt-ncighbor pixels 
can generate a trigger, while signals below the thresh- 
old or pixels not situated closely together cannot con- 
tribute to the decision. This deficiency is particularly 
pronounced in the case of shower images of a light con- 
tent close to the threshold, i.e., in the interesting energy 
region below 60 GeV. For this reason, a new trigger sys- 
tem called sum trigger has been developed to explore the 
energy region down to ~ 25 GeV. Details of the new sum 
trigger will be presented in the Section 2.2 below. 



The signals entering the DAQ system are recorded by 
flash analog-to-digital converters (FADCs) with a sam- 
pling rate of 2 GHz. For each event, 50 FADC slices 
are recorded for all the pixels. The details of the DAQ 
system are described in Goebel et al. (2007). For the pul- 
sar study, the central pixel of the camera was modified 
to record the optical flux from the object under study, 
i.e., to measure the optical pulsations of the Crab pulsar. 
The details of the central pixel system can be found in 
Lucarelli et al. (2008). The telescope tracked the Crab 
position with a typical precision of 0?02. In addition, we 
regularly recorded calibration and pedestal events with a 
frequency of 25 Hz each. Further details on the telescope 
can be found in Baixeras et al. (2004). 

2.2. The Sum Trigger 

As mentioned above, the standard trigger scheme is not 
very efficient below 60 GeV because even at 25 GeV the 
image covers well over four pixels and the signals show 
a wide spread in amplitude. In order to improve the 
trigger efficiency just above threshold, a new technique 
was developed, the so-called sum-trigger method. The 
main feature of the sum trigger is the summation of the 
analog pixel signals from a wider camera area, so-called 
patches, followed by a discrimination of this summed- up 
signal. There are 24 partially overlapping patches in an 
annulus with inner and outer camera radii of <~ 0?25 and 
~ 0?8 respectively. Each patch comprises 18 pixels. The 
threshold level for the summed signal from a patch is an 
amplitude of 27 PhEs. 

This sum-trigger concept has some clear advantages 
compared to the standard trigger. The summation of 
the analog signals allows any pixel signal in the patch to 
contribute to the trigger, even if its amplitude is below 
the pixel threshold of the standard trigger. The con- 
cept, however, has the disadvantages of being quite sen- 
sitive to accidental triggers from afterpulses. Aftcrpulses 
are caused by PhEs hitting the dynodes and sometimes 
releasing ions from adsorbed water or adsorbed gases. 
These ions are back-accelerated by the relatively high 
voltage between the photocathode and the first dynode, 
hit the photocathode and liberate many secondary elec- 
trons. The afterpulsc amplitude spectrum follows ba- 
sically an exponential distribution, which drops signifi- 
cantly slower than the Poissonian night sky light distri- 
bution and completely dominates the rate of signal above 
5-6 PhE. In a single patch, i.e., from the sum of 18 pix- 
els, the rate of afterpulses above 27 PhEs was found to 
be around 20 — 30 kHz, which is far beyond the MAGIC 
DAQ rate limit of 1 kHz. To suppress this undesirable 
background, individual pixel signals were, before sum- 
ming, clipped at an amplitude of <~ 6 PhE, thus making 
the trigger insensitive to large aftcrpulses of individual 
PMTs. The sum-trigger area (the annulus with inner 
and outer radii of ~ 0?25 and <~ 0?8), the patch size (18 
pixels), the threshold level (27 PhE in amplitude after 
sum) and the clipping level (<~ 6 PhE) were optimized 
by measurements and detailed Monte Carlo (MC) stud- 
ies (Rissi 2009). The chosen settings for the sum trigger 
result in a trigger threshold of 25 GeV for a gamma-ray 
source with an index of —2.6 as shown in the top panel of 
Figure 1 (according to the convention in ground-based 7- 
ray astronomy the threshold is defined as the peak of the 
reconstructed differential energy spectrum). The sum 
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Fig. 1. — Comparison of the energy distribution of triggered 
events (top) and the gamma-ray collection area (bottom) for 
MAGIC I between the standard trigger and the sum trigger com- 
puted by MC simulations. For the gamma-ray energy spectrum a 
power law with an index of —2.6 was assumed in the simulations. 

trigger also improved significantly the collection area for 
low energy showers when compared to the area of the 
standard trigger. At 20 GeV the collection area of the 
sum trigger is 10 times larger and at 60 GeV still twice 
as large when compared to that of the standard trigger, 
as shown in the bottom panel of Figure 1. A more de- 
tailed description of the sum trigger is presented in Rissi 
(2009) and Rissi et al. (2009). 

3. OBSERVATION AND DATA PROCESSING 

3.1. Observation 

The first observation of the Crab pulsar with the sum 
trigger started on 2007 October 21 and extended up to 
2009 February 3. In total the Crab pulsar was observed 
for 48 hr in winter 2007/2008 (Aliu et al. 2008) and for 
78 hr in winter 2008/2009. In the first year, all obser- 
vations were restricted to zenith angles below 20° where 
the air mass between the showers and the telescope is 
the lowest possible, i.e., the atmospheric transmission 
for Cherenkov light is highest. In this zenith angle range 
the correlation between energy and observed number of 
PhEs is almost independent of the zenith angle and the 
trigger threshold is nearly constant as a function of the 
zenith angle. In the second year, some of the observations 
were done at zenith angles above 20°. These data are not 
used in the following analysis. In the second campaign, 
five sub-patches malfunctioned. The losses are estimated 
and corrected using MC simulations. 

3.2. Data Processing 

In the calibration process, the conversion factors from 
the FADC counts to the number of PhEs and the relative 



timing offsets of all pixels are computed using the cali- 
bration and pedestal events. The details of the procedure 
can be found in Albert et al. (2008a). After the calibra- 
tion, an image cleaning is performed in order to remove 
pixels which do not contain a useful Cherenkov photon 
signal (e.g., pixels only containing FADC pedestal, NSB 
photons, and afterpulses) . The standard procedure of 
the image cleaning can be found in Aliu et al. (2009). 
Since this study aims for the lowest possible threshold, 
a more sophisticated method of image cleaning was used 
in this analysis. At first the algorithm searches for the 
core pixels of the shower image. The definitions of the 
core pixels are as follows. 

• If two neighboring pixels have more than 4.7 PhE 
each and the time difference is less than 0.8 ns, 
these two pixels are core pixels. 

• If three neighboring pixels have more than 2.7 PhE 
each and the arrival times of all three are within 0.8 
ns, these three pixels are core pixels. 

• If four neighboring pixels have more than 2.0 PhE 
and the arrival times of all four are within 1.5 ns, 
these four pixels are core pixels. 

After the core search, boundary pixels of the image are 
selected. For a pixel to be defined as a boundary pixel, 
the following three conditions must be fulfilled: 

1 . The pixel must be a neighbor to at least one of the 
core pixels. 

2. The pixel must have more than 1.4 PhE. 

3. The time difference to at least one of the neighbor- 
ing core pixels must be less than 1.0 ns. 

The charge and timing information of pixels which are 
neither "core" nor "boundary" is discarded. In this way, 
accidental trigger events (half of the recorded data) are 
efficiently rejected and the contamination of e.g., NSB 
photons to the shower image can be mostly suppressed. 
Further details of this method can be found in Shayduk 
et al. (2005). 

After the image cleaning, the conventional image pa- 
rameters are calculated in the standard way (Hillas 1985; 
Aliu et al. 2009). 

3.3. Data Pre-selection 

Only data taken under stable atmospheric and hard- 
ware conditions were used in the analysis. Selection cri- 
teria include, for instance, the performance of the mirror 
focusing and reflection, the cloud coverage and the sta- 
bility of the event rate after image cleaning. A detailed 
description of the pre-selection criteria can be found in 
Saito (2010). As mentioned before, only data with zenith 
angles below 20° were used in the analysis. After apply- 
ing all criteria 25 hr remained from the winter 2007/2008 
period and 34 hr from the winter 2008/2009 period. 

3.4. Event Selection 

A cleaned image of a 25 GeV 7-ray on average has 
~ 8 pixels, which is barely sufficient to perform a mo- 
ment analysis to obtain the conventional Hillas image 
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Fig. 2. — ALPHA distribution of 7-ray MC events as a function 
of SIZE. Red stars indicate the cut values which maximize the 
Q-factor defined by Equation (1). The black solid line shows the 
STZiJ-dependent ALPHA cut used in the analysis. The green 
dotted line denotes SIZE = 25 PhE, below which data are not 
used for the analysis. 

parameterization of the shower. However, the only ef- 
fective way to separate 7-rays from the background be- 
low 100 GeV exploits the orientation of the images in 
the camera plane. 7-rays from the source have images 
that point with their major axis towards the source lo- 
cation. Background events (mainly hadron showers and 
muon arcs/rings), on the other hand, have images that 
are randomly oriented. The parameter that describes the 
orientation of an image in the camera is called ALPHA, 
which is the angle between the major axis of an image 
and the line connecting the center of gravity (COG) of 
that image with the location of the source in the camera. 
The ALPHA distribution of the 7-ray MC events as a 
function of SLZE is illustrated in Figure 2, where SIZE 
is the total number of PhE in the image. Red stars indi- 
cate the cut values which maximize the so-called quality 
factor Q defined as 



Q 



(1) 



where e 7 and €bg are the survival efficiencies of the 7-ray 
events and hadron background events, respectively. 

In the analysis we used an ALPHA cut depending on 
SIZE, which was derived by fitting a numerical function 
a(\og w {SIZE) + b) c to the best cut values found in the 
individual SIZE bins (stars in Figure 2). e 7 , cbg and 
Q for the used SI ZE- dependent ALPHA cut are shown 
in Figure 3. 

In addition, the so-called spark-like events, which are 
created by the discharge of charge accumulated at the 
glass envelope of some PMTs, are rejected. These events 
are efficiently identified by applying the condition 

4.0 x log 10 {CONC) > 1.5 - \og w {SIZE) (2) 

where, CONC is defined as the sum of the charges in the 
two pixels with the highest content divided by SIZE, 
which indicates how strongly the charge is concentrated 
in a small region. 

3.5. Energy Estimates 
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Fig. 3. — 7-ray survival efficiency (red), hadron background sur- 
vival efficiency (green) and the Q-factor (blue, Equation (1)) as a 
function of SIZE for the used 5i"Z_E-dependent ALPHA cut. 



The energy of each event is estimated by means of 
the Random Forest method. After the training with 
MC 7-ray events, the Random Forest assigns the most 
probable energy to each event by using several image 
parameters in a comprehensive manner. The details 
of the method can be found in Albert et al. (2008b). 
In this analysis SIZE, LENGTH, DIST and the 
zenith angle were used for the Random Forest energy 
estimation. LENGTH is the second moment of the 
charge distribution along the major axis of the shower 
image, while DIST is the distance between the source 
position in the camera and the center of gravity of the 
image. WIDTH (the second moment along the minor 
axis of the image) and CONC , which are normally 
included, were not used because it was found that they 
do not contribute to the energy estimate in the very low 
energy regime. 



3.6. Pulse-phase Calculation 

Each event is marked with a time stamp, which gives 
the time when the event was triggered. The time stamps 
are derived from a GPS controlled Rubidium clock and 
have an absolute accuracy of less than 1 /us. In order 
to compensate the varying propagation times of the 7- 
rays within the solar system, which are mainly due to 
Earth's movement around the Sun, the recorded times 
are transformed to the barycenter of the solar system. 
The barycentric correction was done with the software 
package TEMPO [24]. The rotation frequency vq, its 
time derivative v$ and the barycentric times of the main 
pulse peak of the Crab pulsar are monitored in radio at 
610 MHz by the Jodrell Bank radio telescope (Lyne et al. 
1993) and the values of the parameters are published 
every month 35 . Based on these, the rotational phase of 
each event is computed from the barycentric times with 
the following formula 



Phase = v (t - i ) + -i>o(t - t f 



(3) 



http:/ /www.jb. man. ac.uk/ pulsar/crab. html 
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Fig. 4. — Pulse profile of the Crab pulsar seen in optical wave- 
lengths with 15 hr of observations. The signal was recorded with 
the central pixel of the MAGIC camera in parallel to the 7-ray 
measurements. 

The second and higher derivative terms of this Taylor 
series are negligible on the scale of one month. 

4. PULSATION OF THE CRAB PULSAR ABOVE 
25GEV 

4.1. Optical Pulsation 

In order to assure that the timing system of MAGIC 
and the pulse phase calculation worked properly, the op- 
tical pulsation of the Crab pulsar was checked first. The 
optical pulsation was measured with the central pixel, 
which was modified for this purpose to be sensitive and 
to digitize the light flux variation from the source. Ev- 
ery time a shower event was triggered, the signal of the 
central pixel was recorded by the DAQ for 25 ns. 

The phase distribution (hereafter pulse profile) of the 
central pixel data is shown in Figure 4. Two peaks are 
clearly visible at the expected phases. Phase corre- 
sponds to the main peak position in radio at 610 MHz. 
A delay of ~ 0.01 in phase can be seen with respect to 
the radio main peak position, which is known and con- 
sistent with other measurements (see, e.g., Oosterbroek 
et al. (2008)). 

4.2. Pulsation above 25 GeV 

The pulse profile of the 7-ray events detected with 
MAGIC is shown in Figure 5. Events with SIZE below 
25PhE and above 500 PhE were discarded. Note that 
every event is shown twice (three times for the first bin) 
in order to generate a pulse profile that spans the phase 
region from -45/44 (-1.0227) to 45/44 (1.0227). The 
bin width is 1/22 in phase, which corresponds to about 
1.5 ms. An excess is evident in the profile at the position 
of the main pulse and interpulse of the pulsar. Follow- 
ing the often-used convention (Fierro et al. 1998) of PI 
(main pulse phase from —0.06 to 0.04) and P2 (interpulse 
phases from 0.32 to 0.43), the numbers of excess events 
in PI and P2 are 6200 ± 1400 (4.3 a) and 11300 ± 1500 
(7.4tr), respectively. By summing up PI and P2, the 
excess corresponds to 7.5o\ The background level was 
estimated using the so-called off pulse (OP) phases (0.52 
-0.88, Fierro et al. (1998)). Above 25 GeV, the flux of P2 
is nearly twice that of PI. The width of the main pulse 
is significantly smaller than the conventional PI phase 
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Fig. 5. — Pulse profile of the Crab pulsar obtained with MAGIC. 
The red shaded area indicates the signal phases (PI and P2) while 
the black shaded area indicates the background control phases (OP 
phases). 
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Fig. 6. — Pulse profile of the Crab pulsar when the ALPHA cut 
is inverted. As expected, the signals seen in Figure 5 disappear. 

interval. By defining the signal phase of the main pulse 
to be —1/44 to 1/44 in phase (the bin with the largest 
number of events), the excess is 6400 ±970 corresponding 
to 7.0cr. No significant emission between the main pulse 
and the interpulse was detected. A detailed study of the 
pulse profile is given in Section 7. 

In order to verify the soundness of the signal, further 
tests were made. Firstly, the phase distribution of the 
events which are rejected by the ALPHA cut is exam- 
ined. The results are shown in Figure 6. As expected, the 
distribution is consistent with statistical fluctuations of 
the background without any signal. Also, the growth of 
the number of excess events as a function of the number 
of background events is checked. As one can see in Fig- 
ure 7, the excess grows linearly, assuring that the signal 
is constantly detected. 

4.3. Variability Study 

The linear growth of the excess shown in Figure 7 im- 
plies a constant flux of the pulsed signal. Nevertheless, 
we also applied the % 2 method to test for a possible yearly 
variability. The number of excess events as a function of 
SIZE is compared between the two years in Figure 8. 
The difference of observation time and the effect of the 
malfunctioning of the sum-trigger sub-patches are cor- 
rected for the second year data. Using MC simulations, 
the sub-patch malfunction effects on the acceptance were 
estimated to be 21%, 17%, 11% and 7% for the SIZE 
ranges of 25-50, 50-100, 100-200 and 200 - 400 PhE, re- 
spectively. The x 2 values for the comparison of the two 
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Fig. 7. — Numbers of excess events as a function of the number of 
background events (events in OP phases) for PI (red), P2 (blue), 
and the sum of the two (black). They grow linearly, implying 
constant detection of the signal. 

years are 1.0 and 3.1 for PI and P2, respectively, for 4 
degrees of freedom. No significant yearly variability was 
found in the flux. 

We also studied a possible variability of the pulse pro- 
file. Figure 9 shows the pulse profiles for the first (win- 
ter 2007/20088) and the second (winter 2008/2009) year. 
The two profiles are compared with each other with a 
X 2 test from phase -0.0682 to 0.432 (-3/44 to 19/44), 
which is roughly from the beginning of PI to the end of 
P2. The obtained x 2 is 5.0 for 10 degrees of freedom. 
Therefore, no significant yearly variability of the pulse 
profile was found. 
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Fig. 8. — The comparison of the SIZE distributions between 
the winter 2007/2008 and the winter 2008/2009. The difference in 
observation time and the hardware performance are corrected by 
scaling the data of winter 2008/2009. 



5. EVALUATION OF THE EXPONENTIAL 
CUTOFF SPECTRUM SUGGESTED FROM THE 
FERMI-LAT DATA 

It is very important for the verification of the stan- 
dard OG and SG models to check if the energy spectrum 
follows an exponential cutoff. All energy spectra mea- 
sured by Fermi-LAT up to a few tens of GeV are indeed 




Fig. 9. — Pulse profiles of the Crab pulsar for the winter 
2007/2008 (top) and the winter 2008/2009 (bottom). The two 
profiles are statistically consistent. 

consistent with the OG/SG model. In this section, we 
evaluate the exponential cutoff hypothesis based on the 
measurements performed by Fermi-L AT and MAGIC. 

5.1. Analysis of Public Fermi-L AT Data 

Although the Fermi-LAT Collaboration published 
their results of the Crab pulsar observations (Abdo et al. 
2010a), we performed a customized analysis of the spe- 
cific phase intervals in order to properly compare the 
Fermi-LAT and MAGIC data. One year of FermirLAT 
data taken from 2008 August 4 to 2009 August 3 are 
analyzed. Events with an energy between 100 MeV and 
300 GeV and with an arrival direction within a radius 
of 20 degrees around the Crab pulsar were downloaded 
from the Fermi Science Support Center. 36 Only events 
with the highest probability of being photons, those in 
the diffuse class, were used in this analysis. Events with 
imperfect spacecraft information and events taken when 
the satellite was in the South Atlantic Anomaly were 
rejected. In addition, a cut on the maximum zenith 
angle (< 105°) was applied to reduce the contamina- 
tion from the Earth-albedo 7-rays, which are produced 
by cosmic rays interacting with the upper atmosphere. 
The pulse-phase assignment to each event was carried 
out by the Fermi-LAT analysis tool with the monthly 
ephemeris information from Jodrell Bank. 37 For the de- 
tector response function, "P6_V3_Diffuse" is used, while 
"isotropic_iem_v02.txt" and "gll_iem_v02.fit" are used for 
the extragalactic and galactic diffuse emission models. 
In order to estimate the contribution of the Crab Neb- 
ula, the pulse-phase interval between 0.52 and 0.87 is 
used. 38 . The contamination from nearby bright sources 
such as Geminga and IC 443 are also taken into account 

36 http:/ /fermi. gsfc.nasa.gov/ssc/data/ 

37 http://www.jb.man.ac.uk/ pulsar/crab. html 

38 This phase range is not identical to the OP phases (0.52 — 0.88; 
Fierro et al. (1998)), which was used for the MAGIC data to esti- 
mate the background level. We adopt this range to be consistent 
with Abdo et al. (2010a) 
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Fig. 10. — Top: the pulse profile of the Crab pulsar above 
100 MeV produced with the Fermi- LAT data. Bottom: the en- 
ergy spectrum of the Crab pulsar (total pulse). The black line and 
dots are obtained from the public Fermi-LAT data, while the green 
line is the spectrum reported in Abdo et al. (2010a). 

in the calculation of the spectrum. The unbinned like- 
lihood spectral analysis assuming a power-law spectrum 
with an exponential cutoff 



d 3 F(E) 
dEdAdt 



F^E/lGeYy^eM-E/Ec) (4) 



gives Fx = (2.32 ± 0.05 stat ) x 1CT 10 cm^s^MeV- 1 , 
ri = 1.99 ± 0.02 stat and E c = 6.1 ± 0.5 stat GeV as 
best fit parameters, consistent with the values in Abdo 
et al. (2010a). Results for the total pulse are shown by 
the solid black line in Figure 10, together with the pulse 
profile above 100 MeV. The green curve in the figure is 
the spectrum given in Abdo et al. (2010a). The points are 
obtained by applying the same likelihood method in the 
limited energy intervals, assuming a power-law spectrum 
within each interval. 

The same analysis was applied to PI, P2 and the sum 
of the two (PI + P2). The results are shown in Figure 
12 and the best-fit parameters are summarized in Table 
1. 

5.2. Statistical Evaluation of the Difference between 
the Extrapolated Fermi-LAT Spectrum and the 
MAGIC Data 

For MAGIC, the energy resolution below 50GeV is 
about 40%. In addition, near the trigger threshold, the 
energy is overestimated because only events which de- 
posit more Cherenkov photons onto the MAGIC mirrors 
are selectively triggered. These effects can be corrected 
with MC simulations, but the correction introduces ad- 
ditional systematic uncertainties. In order to minimize 
these uncertainties, we adopt the following method. As- 
suming that the exponential cutoff spectrum determined 
by the Fermi-LAT below a few tens of GeV is valid up to 
2TeV, we calculate the expected SIZE distribution in 



a, 9000 
c 

£8000 
7000 
6000 
5000 
4000 
3000 
2000 
1000 




«4000 
c 

£3500 



3000 
2500 
2000 
1500 
1000 
500 


,6000 



LU 

*5000 



4000 



3000 



2000 



1000 



P1 +P2 



Expected 

Expected (E:+23%; 
Measured 



10^ 



SIZE [PhE] 



PI 



Expected 

Expected (E:+23%; 
Measured 



10^ 



SIZE [PhE] 



P2 




Expected 

Expected (E:+23%; 
Measured 



10^ 



SIZE [PhE] 



Fig. 11. — Comparison of the SIZE distribution between the ex- 
pectations and the measurements. From the top, PI + P2, PI and 
P2 are shown. Red lines show the measured distributions while 
the solid black lines indicate the expected distributions computed 
with MC simulations assuming the exponential cutoff spectra de- 
termined by Fermi-LAT (see Table 1 and Figure 12). Dotted lines 
are the expectation in the case where the energy scale of MAGIC 
is 23% higher than that of Fermi-LAT due to systematic uncer- 
tainties of both instruments. 
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the MAGIC data. Taking into account the trigger thresh- 
old of 27 PhE, the events with SIZE below 30 PhE are 
not used to avoid a possible mismatch between MC and 
real data near the threshold. Then, we compute the % 2 
value between the expected and measured distributions. 
Statistical errors of E c measured by F ermi-LAT are taken 
into account as an error of the expected distribution. 

The results are shown in Figure 11. The x 2 values are 
54.2 , 15.8 and 42.3 for PI + P2, PI and P2, respectively. 
The number of degrees of freedom is 3. The exclusion 
probabilities correspond to 6.7a, 3. Oct and 5.8ct. It should 
be noted that the possible energy scale shift between the 
two instruments is not taken into account here. 

The systematic uncertainty in the energy scale of the 
Fermi-LAT is estimated to be less than 7% above 1 GeV 
(Abdo et al. 2009) while that of MAGIC is estimated 
to be 16% (Albert et al. 2008c). We performed the same 
statistical test with an increased Fermi-LAT energy scale 
of 23%, corresponding to the linear sum of the system- 
atic errors in both instruments. The results are shown as 
dotted lines in Figure 11. Though the discrepancies be- 
come smaller, the x 2 values are 42.3, 12.6 and 30.0 with 
the number of degree of freedom 3 for PI + P2, PI and 
P2, respectively. The exclusion probabilities correspond 
to 5.8ct, 2.5ct and 4.7ct. Even with the systematic uncer- 
tainties taken into account, the inconsistency between 
the extrapolated Fermi-LAT spectrum and the observa- 
tions by MAGIC is significant. 

6. ENERGY SPECTRA 

In the previous section, it was shown that the ex- 
trapolation of the Fermi-LAT measured spectra under 
the exponential cutoff assumption results in significant 
differences with the MAGIC data above 25 GeV. Here 
we present the energy spectrum between 25 GeV and 
100 GeV based on the MAGIC measurements, which 
are the first flux measurements in this energy region 
and complement the Fermi-LAT and VERITAS measure- 
ments (Aliu et al. 2011). 

6.1. Spectra of PI, P2 and PI + P2 

The energy spectra of the Crab pulsar were computed 
based on the detected excess events found in PI and P2, 
using the standard MAGIC software. The energy reso- 
lution and the trigger bias effect were corrected by an 
unfolding procedure which includes the Tikhonov regu- 
larization method (Tikhonov & Goncharsky 1987). The 
results are shown in Figure 12. The combined spectrum 
of PI + P2 is consistent with a power law, which can be 
described using the following formula: 



Spectrum: P1 + P2 
^ 10 



d 3 F(F) 
dEdAdt 



= F 30 (F/30 GeV)" 



F 30 = (14.9±2.9 stat ±9.6 syst )xl0- 



(5) 

^TcV- 1 and 

T 2 = 3.4 ± 0.5 s tat ± 0.3 sys t were obtained as best- 
fit parameters. Also, the spectrum seems to connect 
smoothly to the VERITAS measurements above 100 GeV 
(Aliu et al. 2011). 

The individual spectra of PI and P2 were also calcu- 
lated using the same data set. Results are shown as well 
in Figure 12. The best- fit parameters are summarized in 
Table 1. 
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Fig. 12. — Energy spectra of the Crab pulsar for PI + P2, PI, and 
P2 from the top. The black solid lines and dots are obtained from 
the public Fermi-LAT data, while red points denote the MAGIC 
measurements. The pink line and a butterfly shape indicate the 
power-law fit to the MAGIC data and its statistical uncertainty. 
The gray shade with a black line and open diamonds denote the 
VERITAS measurements (Aliu et al. 2011). The dotted lines show 
the results of the combined fit above 5 GeV (see Section 6.2). 

Figure 12 clearly shows the deviation of the MAGIC 
spectra with respect to the extrapolation of the expo- 
nential cutoff spectra determined by Fermi-LAT, which 
is consistent with our statistical analysis in the previous 
section. 

6.2. Combined Fit above 5 GeV 

To get a better estimate of the power-law index for the 
higher energies, Fermi-LAT data points above 5 GeV and 
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TABLE 1 

The best fit parameters of the spectra for different phase intervals. 





Fermi-LAT a 


MAGIC" 


Phase 


Ft [10 _1 "cm _ ^s _1 MoV _1 ] 


Ti 


E c 


F 3 o [lO^cm^s^TcV- 1 ] 


r 2 


Total 


2.32 ± 0.05 stat 


1.99 ± 0.02 stat 


6.1 ± 0.5 stat 






PI + P2 


1.94 ± 0.05 stat 


1.98±0.02 stat 


4.5±0.3 stat 


14.9±2.9 stat ±9.6 SBSt 


3.4±0.5 stat ±0.3 s „ st 


PI 


1.29±0.04 stat 


1.99±0.02 stat 


3.7±0.3 stat 


4.5±2.3 3tat ±2.e syst 


3.1 ± 1.0 stat ± 0.3 s „ st 


P2 


0.67 ± 0.02 stat 


1.95±0.03 3tat 


5.9±0.7 stat 


10.0 ± 1.9 stat ±6.7 SBSt 


3.4±0.5 stat ±0.3 s „ st 



a Obtained by fitting Equation (4) to Fermi-LAT data. 
b Obtained by fitting Equation (5) to MAGIC data. 



MAGIC data points are combined and fitted by a power 
law. 



d 3 F(E) 
dEdAdt 



= Fio(£/10GeV)- 



(6) 

It should be mentioned that Fermi-LAT points are ob- 
tained using the likelihood analysis for each energy in- 
terval assuming a power law, while for the fit each point 
was assumed to follow Gaussian statistics with the stan- 
dard deviation being the error obtained by the likelihood 
results. Though this does not statistically correspond to 
the exact likelihood function of the problem, it is a good 
and appropriate approximation. The results are shown 
in Table 2. 



7. A CLOSER LOOK AT THE PULSE PROFILES 

7.1. Peak Phase and Pulse Width 

We examine the peak phase and the pulse width in the 
MAGIC energy range assuming a pulse shape a priori and 
fitting it to the measured profile. The used functions are 
a Gaussian 



F G (,) = F exp(-^-^) 



and a Lorcntzian 



F L {x) =F 1 + 



(7) 



(8) 



where fi corresponds to the peak phase, while a can be 
translated into a pulse width. The full width at half- 
maximum (FWHM) is equal to 2.35 x a and 2 x a for the 
Gaussian and the Lorentzian, respectively. In order to 
study the asymmetry of the pulses, asymmetric Gaussian 
and Lorentzian functions are also assumed with a being 
different below and above x = /i. 

The results are shown in Table 3. The peak phase 
of the main pulse is compatible with 0.0 (defined by 
the radio peak) for all fitted parametcrizations. The 
Fermi-LAT Collaboration reported that the pulse shape 
is well modeled by an asymmetric Lorentzian and the 
peak phase above 100 MeV is -0.0085 ± 0.0005 (Abdo 
et al. 2010a). The MAGIC result under the asymmet- 
ric Lorentzian assumption is consistent with it as well. 
The peak phase of the interpulse depends on the assump- 
tion of the shape, while it is approximately 0.39 ± 0.01, 
which is also consistent with the value above 100 MeV 



TABLE 2 

The power law fit (Equation (6)) combining the Fermi-LAT 

DATA ABOVE 5 GeV AND THE MAGIC DATA. 



Phase 


Fio [10"' cm- x s- I TcV- 1 ] 


r 


X 2 /n.d.f 


PI + P2 
PI 
P2 


3.0 ± 0.2 
1.2 ± 0.2 
1.7 ± 0.2 


3.0 ± 0.1 

3.1 ± 0.2 
2.9 ± 0.1 


8.1/4 
1.6/4 
8.3/4 



(0.398 ±0.003, Abdo et al. (2010a)). The FWHM of the 
main pulse is approximately 0.03 ± 0.01 independently 
of the assumed shape, while that of the interpulse is 
0.07 ± 0.01 and 0.05 ± 0.01 for the assumption of the 
Gaussian shaped and Lorentzian shaped pulse, respec- 
tively. The main peak is narrower than the interpulse. 
The asymmetric assumptions imply that, for the main 
pulse, the rising edge is steeper than the falling edge, 
while the opposite is true for the interpulse. In Figure 14, 
the half- widths a for the rising and falling edges of both 
the main pulse and the interpulse are compared with the 
values reported in Abdo et al. (2010a). The rising half of 
the main pulse become narrower as the energy increases. 
Though the uncertainty is larger, a similar tendency is 
also visible in the rising and falling halves of the inter- 
pulse, while no such energy dependence is visible in the 
falling half of the main pulse. 

7.2. Other Emission Components 

The AGILE Collaboration reported a possible third 
peak at phase between 0.65 and 0.8 above 100 MeV with 
the significance of 3.7cr (Pellizzoni et al. 2009). A hint 
of a third peak is also seen in the Fermi-LAT data at 
phase ^0.74 only above 10 GeV with the significance of 
2.3cr (Abdo et al. 2010a). They coincide with the radio 
peak observed between 4.7 and 8.4 GHz. In the MAGIC 
data above 25 GeV, a similar peak is seen at phase ~ 

0. 73. Defining the signal phases as 0.72 — 0.75 and the 
background control phases as the OP phases (0.52 — 0.88) 
excluding the signal phases, 1600±700 excess events were 
found, corresponding to 2.2a\ This pre-trial significance 
is too low to claim a detection and it is within the range 
of expected fluctuation of the background. 

An emission between the main pulse and the interpulse, 

1. e., a so-called bridge emission is seen in some energy 
bands. With the MAGIC data above 25 GeV, it is not 
visible though the statistical uncertainty is large. Defin- 
ing the signal phases as 0.04 — 0.32, i.e., from the end of 
PI to the beginning of P2, and using the OP phases as 
the background estimate, 3200 ±2800 excess events were 
found, corresponding to 1.1c The upper limit on the 
number of excess events with the 95% confidence level is 
8800, which corresponds to half of the flux of PI + P2. 
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Fig. 13. — Pulse profile of the MAGIC data fitted by asymmetric Lorentzians (red lines) and asymmetric Gaussians (black lines). 
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Fig. 14. — Energy dependence of the rising and falling half-width for the main pulse (left) and the interpulse (right) assuming an 
asymmetric Lorentzian as the pulse shape. Points below 10 GeV are reproduced from Abdo et al. (2010a). 



TABLE 3 

Fit results using different assumptions on the shape of the main pulse and the interpulse. 





Main pulse 


Assumed Shape 








FWHM 


X^/n.d.f 


Gaussian 


0.7 ±2.6 


14.2 ± 2.4 




33.3 ±5.7 


13.6/17 


Asym. Gauss. 


-2.2 ±6.6 


11.8 ±5.4 


16.9 ±7.0 


33.7 ± 10 


13.3/16 


Lorentzian 


0.4 ±2.7 


14.4 ±3.8 




28.8 ±7.6 


13.4/17 


Asym. Lorentz. 


-8.9 ±6.5 


5.9 ±6.3 


25.4 ±8.0 


31.3 ±10 


11.2/16 




interpulse 


Assumed Shape 


/" 






FWHM 




Gaussian 


377.2 ±5.7 


32.4 ±4.8 




76 ± 11 


38.5/37 


Asym. Gauss. 


391.8 ± 10 


42.1 ± 9.1 


18.8 ±8.2 


72 ± 14 


36.6/36 


Lorentzian 


384.1 ±5.2 


26.2 ± 7.2 




52 ± 14 


41.6/37 


Asym. Lorentz. 


392.9 ±8.7 


37.2 ± 12.2 


15.6 ±8.2 


53 ± 15 


39.5/36 



Units of all the parameters are 10 3 in phase 
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8. DISCUSSION AND CONCLUSIONS 

8.1. Summary of Findings 

The findings of this study can be summarized as fol- 
lows. 

1. 59 hr of MAGIC observations of the Crab pulsar 
during the winters 2007/2008 and 2008/2009 re- 
sulted in the detection of 6200 ± 1400 and 11300 ± 
1500 excess events from PI and P2, respectively. 
The flux of P 2 is a factor ~ 2 larger than that of 
PI above 25GcV. 

2. No yearly variability in the pulse profile or in the 
flux was found. 

3. The flux measured with MAGIC is significantly 
higher than the extrapolation of the exponential 
cutoff spectrum determined by Fermi-LAT. 

4. The energy spectra extend up to at least 100 GeV 
and can be described by a power law between 
25 GeV and 100 GeV. The power-law indices of 
PI, P2 and PI + P2 are -3.4 ± 0.5 sta t ± 0.3 syst , 
-3.1±1.0 stat ±0.3 syst and -3.4±0.5 stot ±0.3 S!/s t, 
respectively. The sensitivity of MAGIC I above 
100 GeV is not sufficient to clarify if the spectrum 
continues with a power law or drops more rapidly. 
However, our spectrum of PI + P2 and the VER- 
ITAS measured spectrum above 100 GeV seem to 
be a good extrapolation of each other. 

5. Assuming an asymmetric Lorcntzian for the pulse 
shape, the peak positions of the main pulse and the 
interpulse are estimated to be — 0.0089±0.0065 and 
0.3929 ± 0.0087 in phase, while the FWHMs arc 
0.031 ± 0.010 and 0.053 ± 0.015. Compared with 
the Fermi-LAT measurements, the pulse widths are 
narrower in the MAGIC energy regime. 

6. The bridge emission between PI and P2 is weak. 
With the current sensitivity no signal was found. 
A potential third peak with a pre-trial significance 
of 2.2cr is seen at a similar position as in the AG- 
ILE data above 100 MeV and the Fermi-LAT data 
above 10 GeV, but it is consistent with the back- 
ground fluctuation. 

The spectrum of the Crab pulsar does not follow an 
exponential cutoff but, after the break, it continues as 
a power law. This is inconsistent with the OG and SG 
models in their simplest version, where it is assumed that 
the emission above <~ 1 GeV comes only from curvature 
radiation, leading to an exponential cutoff in the spec- 
trum. 

A theoretical interpretation of this deviation from the 
exponential cutoff is discussed in the next section. 

8.2. Theoretical Interpretation of the Spectrum 

In a pulsar magnetosphere, high-energy photons are 
emitted by the electrons and positrons that are accel- 
erated by the magnetic-field-aligned electric field, E\\. 
To derive E\\ , we must solve the inhomogeneous part of 
the Maxwell equations (Fawley et al. 1977; Scharlemann 
et al. 1978; Arons & Scharlemann 1979) 

V • E\\ = 4n(p - p GJ ) (9) 



where p denotes the real charge density, and pqj = 
Q, ■ B/(2nc) the Goldreich— Julian charge density (Gol- 
dreich & Julian 1969; Mestel 1971); the angular- velocity 
vector ft points in the direction of the NS spin axis with 
magnitude |fi| = 2tt/P, B refers to the magnetic field, 
and c is the speed of light. If p coincides with pqj m the 
entire magnetosphere, -Ey vanishes everywhere. However, 
if p deviates from pcj m some region, it is inevitable for 
a non-vanishing E\\ to arise around the region and the 
particle accelerator, or the so-called the "gap" , appears. 

To predict the absolute luminosity of the gap, as well 
as any phase-averaged and phase-resolved emission prop- 
erties, we must constrain = |2?|||, p, and the gap 
geometry in the three-dimensional magnetosphere. We 
solve the Poisson equation (9) together with the Boltz- 
mann equation for electrons and positrons (e ± 's) and 
with the radiative transfer equation between 0.005 eV 
and 10 TeV. The position of the gap is solved within the 
free-boundary framework so that may vanish on the 
boundaries, and turned out to distribute in the higher al- 
titudes as a quantitative extension of previous OG mod- 
els, which assume a vacuum (i.e., p = 0) in the gap. 
In the present paper, we propose a new, non-vacuum 
(i.e., p 7^ 0) OG model, solving the distribution of p self- 
consistently from the e~-e + pair creation at each point 
(Beskin et al. 1992; Hirotani & Okamoto 1998; Hirotani 

6 Shibata 1999a,b; Hirotani et al. 1999; Takata et al. 
2006). 

The created e ±5 s are polarized and accelerated by E\\ 
in the gap to attain high Lorentz factors, ~ 10 7 ' 5 . Such 
ultra-relativistic e + 's and e _ 's emit primary 7-rays via 
synchro-curvature and inverse Compton (IC) processes. 
For the IC process, the target photons are emitted from 
the cooling NS surface, from the heated PC surface, and 
from the magnetosphere in which pairs are created. The 
primary 7-rays that are emitted by the (inwardly accel- 
erated) e~'s efficiently (nearly head-on) collide with the 
surface X-rays to materialize as the primary e s within 
the gap and as the secondary e ± 's outside the gap. The 
secondary e ±5 s efficiently lose their energy via the syn- 
chrotron process in the inner magnetosphere and cascade 
into tertiary and higher-generation e ±! s via two-photon 
and one-photon (i.e., magnetic) pair-creation processes. 
The primary 7-rays that are emitted by the (outwardly 
accelerated) e + 's via the curvature process collide with 
the magnetospheric X-rays to materialize as secondary 
e s with initial Lorentz factors 7 ~ 10 3 ' 5 outside the 
gap, while those emitted via the IC process collide with 
the magnetospheric IR— UV photons to materialize with 

7 <~ 10 7 . The secondary pairs with 7 ~ 10 3 5 emit syn- 
chrotron emission below 10 MeV and synchrotron-self- 
Compton (SSC) emission between 10 MeV and a few 
GeV (by up scattering the magnetospheric X-ray pho- 
tons). The secondary pairs with 7 ~ 10 7 , on the other 
hand, emit synchrotron emission below 10 GeV and SSC 
emission between 10 GeV and 1 TeV (by up-scattering 
the magnetospheric IR— UV photons). Such secondary 
SSC photons (between 10 GeV and 1 TeV) are efficiently 
absorbed colliding with the magnetospheric IR— UV pho- 
tons to materialize as tertiary pairs with 10 4 < 7 < 10 6 . 
In the high-energy end (i.e., 10 5 < 7 < 10 6 ), the tertiary 
pairs up-scatter magnetospheric IR— UV photons into 
the energy range between 1 GeV and 300 GeV There- 
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Fig. 15. — Phase-averaged spectrum of the pulsed emission from 
the Crab pulsar predicted by the self-consistent outer Gap model. 
The thin solid line (labeled "1") represents the flux of the photons 
emitted by the primary positrons accelerated in the gap, while the 
thin dashed one (labeled "2") and the thin dotted one (labeled 
"3") by the secondary and tertiary pairs, respectively, created out- 
side the gap. The thick green solid line includes magnetospheric 
absorption and subsequent reprocesses, and hence represents the 
flux to be observed. Interstellar absorption is not considered. The 
filled circles (LECS), open circles (MECS), filled triangles (PDS) 
denote the BeppoSAX observations, while the open triangles the 
Gamma-ray Imaging Spectrometer (GRIS). Inverse filled triangles 
(OSSE), open diamonds (COMPTEL) denote CGRO observations. 
The filled squares denote the Fermi-LAT observations, while the 
red filled diamonds the MAGIC observations (this work) . The but- 
terfly shape above 100 GeV indicates the VERITAS observations 
(Aliu et al. 2011). The ordinate is in (MeV) 2 s" 1 cm' 2 MeV -1 
unit. Data points are from Kuiper et al. (2001) and Abdo et al. 
(2010a) 

fore, if we focus the emission component above 10 GeV, 
the photon are emitted by the secondary and tertiary 
pairs when they up-scatter the magnetospheric IR— UV 
photons (see also Lyutikov et al. (2011) for an analytical 
discussion of this process using the multiplicity factor of 
higher- generation pairs). Note that the primary IC pho- 
tons are emitted from relatively inner part of the outer 
magnetosphere and hence totally absorbed by the magne- 
tospheric IR-UV photon field, of which specific intensity 
is self-consistently solved together with En and the par- 
ticle distribution functions at each point. Note also that 
the surface X-ray field little affects the pair creation or 
the IC process in the outer magnetosphere of the Crab 
pulsar. For the details of this self-consistent approach, 
see Hirotani (2006). 

In Figure 15, we present the solved spectral energy 
distribution of the total pulse component when the mag- 
netic axis is inclined 60° with respect to the rotational 
axis and when the observer's viewing angle is 75°. The 
temperature of the cooling NS emission is assumed to be 
70 eV. The distance to the pulsar is assumed to be 2 kpc. 
The thick, green solid line represents the spectrum to be 
observed (i.e., with magnetospheric absorption), while 
the thin black solid line indicates the photons emitted 
by the primary positrons (with Lorentz factor ~ 10 7 ' 5 in- 
side the gap), and the thin dashed line does those emitted 
by the secondary pairs (with Lorentz factor ~ 10 3 ' 5 or 
10 7 outside the gap). Interstellar absorption is not taken 
into account. The primary, un-absorbed IC component 
becomes prominent above 40 GeV; however, most of such 
photons are absorbed by two-photon pair production and 



reprocessed in lower energies as the secondary SSC com- 
ponent. In the secondary emission component, which is 
depicted by the dashed curve, there is a transition of the 
dominant component: the synchrotron component domi- 
nates the IC component (due to the SSC process) below a 
few GeV, whereas the latter dominates the former above 
this energy. 

It should be noted that the pulsed emission between 
25 GeV and 180 GeV is dominated by the SSC compo- 
nent emitted by the secondary and tertiary e^'s. Since 
the higher-generation components, which are denoted by 
the dotted curves, arc emitted from the higher altitudes 
(near the light cylinder), they are less efficiently ab- 
sorbed, thereby appearing as pulsed flux above 20 GeV. 
The resultant spectrum (green, thick solid line) exhibits a 
power-law-like shape above 20 GeV, rather than an expo- 
nential cutoff. We, therefore, interpret that the detected 
7-rays above 25 GeV are mainly emitted via the SSC 
process when the secondary and tertiary pairs up scatter 
the magnetospheric synchrotron IR— UV photons. 

Although the present theoretical result is obtained by 
simultaneously solving the set of Maxwell and Boltz- 
mann equations under appropriate boundary conditions, 
it does not rule out other possibilities such as the syn- 
chrotron emission (e.g., Chkheidze et al. (2011)) or the IC 
emission (e.g., Bogovalov & Aharonian (2000)) from the 
wind zone by ultra-relativistic particles, which may be 
accelerated by MHD interactions or by magnetic recon- 
nection, for instance. However, the fuller study of other 
theoretical models lies outside the scope of the present 
paper. 

8.3. Outlook 

For further studies of the pulsar emission mechanisms, 
observations with a higher sensitivity are essential. The 
stereoscopic system comprising the two MAGIC tele- 
scopes has a 50 GeV threshold and its sensitivity above 
100 GeV is nearly three times higher than that of MAGIC 
I (Aleksic et al. 2011a). The results of stereoscopic ob- 
servations are given elsewhere (Aleksic et al. 2011b). 
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